Orbital Dynamics: The Origin of Anomalous Magnon Softening in Ferromagnetic 

Manganites 
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We study the renormalization of magnons by charge and coupled orbital-lattice fluctuations in 
colossal magnetoresistance compounds. The model considered is an orbitally degenerate double- 
exchange system coupled to Jahn- Teller active phonons. The modulation of ferromagnetic bonds 
by low-energy orbital fluctuations is identified as the main origin of the unusual softening of the 
zone-boundary magnons observed experimentally in manganites. 



PACS number(s): 75.30.Ds, 75.30.Et 

Recently, distinct new features in the spin dynam- 
ics of the ferromagnetic manganese oxide compound 
Pr .63Sr .37MnO3 have been reported Q. Striking de- 
viations from the predictions of the canonical double- 
exchange (DE) theory were observed. In particular, 
an anomalous softening of magnons at the zone boundary 
was found even well below the Curie temperature Tc- It 
is noted that these anomalies are closely related to the 
reduced values of Tc'. For higher- Tc compounds no con- 
siderable deviations from a simple cosine dispersion is 
observed M. 

These experimental findings seem to be of high im- 
portance. They in fact indicate that some very specific 
features of magnetism in colossal magnetoresistance man- 
ganites have still to be identified. In this paper we pro- 
pose a mechanism which might explain the above ex- 
perimental observations. Our basic idea is the follow- 
ing: The strength of ferromagnetic interaction in a given 
bond strongly depends on which orbital is occupied by 
an e g electron (Fig. [I]). Suppose that orbitals and Jahn- 
Teller (JT) distortions are disordered in the ferromag- 
netic phase. It is then evident that temporal fluctuations 
of orbitals affect the short-wavelength magnons through 
a strong modulation of exchange bonds. Quantitatively, 
the effect is expected to be controlled by the character- 
istic time scale of orbital fluctuations: A slowing-down 
of the dynamics of the coupled system of orbitals and 
JT phonons should lead to a stronger magnon renor- 
malization. In other words, the observed zone-boundary 
magnon softening is interpreted in this picture as a pre- 
cursor effect in the proximity of static orbital-lattice or- 
dering. 

More specifically, we calculate the dispersion of one- 
magnon excitation at zero temperature. First, we map 
the ferromagnetic Kondo-lattice model onto a Hamilto- 
nian of interacting magnons and spinless fermions. Yet 
the fermions carry orbital quantum number, the fluctua- 
tions of which the magnon can scatter on. Second, we 
show that as the orbitals are strongly coupled to the 
lattice there is an indirect coupling of magnons to JT 
phonons via the orbital sector. Put another way, the 



orbital-lattice coupling produces a low-energy component 
in the orbital fluctuation spectrum at phonon frequen- 
cies. We then calculate the self-energy corrections to the 
magnon dispersion perturbatively (employing a 1/S ex- 
pansion). Charge fluctuations are found to produce a 
moderate softening of magnons throughout the Brillouin 
zone. The effect of orbital and lattice fluctuations on 
the magnon dispersion is more pronounced and becomes 
dramatic as static order in the orbital-lattice sector is 
approached. 

We start with a model describing the ferromagnetic 
coupling of double-degenerate e g -band electrons to an 
otherwise noninteracting system of localized spins St on 
a cubic lattice: 

H = ~ E *f ( e Le JCT/3 + H.c.) - J H S it s ie 
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+ ^2 Un^ a n. lla + ^2 {u' - JhP) ri ia rnp (1) 

with 7 S {x,y, z} and P = (s ia Sip + |). Hereafter, the 
indices a/f3 and er stand for orbital and spin quantum 
numbers of e g electrons, respectively; summation over 
repeated indices is implied and double counting is ex- 
cluded from the primed sum. The spin operator Si a acts 
on orbital a and Si e = J2a s ia- The last two terms in 
Eq. (Q) describe intra- (inter-) orbital Coulomb interac- 
tion U (U 1 ) and Hund's coupling between e g electrons 
in doubly occupied states. Correlations are assumed 
to be strong: U,U',(U' — Jh) S> t. The important 
point is the peculiar orbital and bond dependencies of 
the electron hopping matrix elements. In orbital basis 
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We work in the limit Jh — > oo. Then the combined ac- 
tion of the on-site Hund's "ferromagnetism" and the spin- 
diagonal nature of electron hopping results in a global 
ferromagnetic ground state. Two different contributions 
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FIG. 1. The e 9 -electron transfer amplitude which controls 
the double-exchange interaction Jde strongly depends on or- 
bital orientation. 



to the spin stiffness should be noticed: Conventional dou- 
ble exchange due to metallic charge motion and superex- 
change which accounts for the kinetic energy gain by vir- 
tual hoppings of e g electrons. We now derive the effective 
Hamiltonian describing the spin excitations at T <C Tq- 
Double exchange. — At Jh — > oo, St and the e g spin 
s e are not independent anymore: Si e — ni e St/2St- They 
form a total on-site spin Si with spin value St + \ if an 
e g electron is present, rii e — 1, or simply St otherwise. 
The unification of band and local spin subspaces implies 
spin-charge separation. That is, the spin component bi a 
of the e g electron Ci aa = bi a Ci a is decoupled and "ab- 
sorbed" by the total spin, and we are left with a charged 
fermion Ci a carrying the orbital index. At T <C Tc, the 
Schwinger boson bti simply becomes a part of the magnon 
operator bi, namely bii — bi/(2S) 1 ^ 2 , while bi\ is almost 
condensed, b^ — (n ie — ri^) 1 / 2 sa 1 — b\bi/AS. Here- 
after, S = St + \. Further, we keep the relevant leading 
1/S terms only and assume small doping x. The kinetic 
energy then reads 



(ij)-t 



(2) 



The first term here describes the fermionic motion in a 
ferromagnetic background (yet this motion is strongly 
correlated in the presence of orbital disorder while 
the second term controls the spin dynamics and spin- 
fcrmion interaction. The latter term would come out 
from a spin-wave expansion of an effective Heisenberg 
model J DE (S l S J ) with J DE - {2S 2 )-Hf (c\ a c ]f3 ) if one 
considers the fermionic sector on average as in mean- 
field treatments of the DE model. However, the bond 
variable c\ a Cjp is a fluctuating complex quantity and the 
spin structure of Eq. (|^) is therefore not of Heisenberg 
form. Explicitly (T < T c ): 
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(3) 



In the classical limit for spins, one recovers from Eq. 
([|) the effective fermionic model with phase-dependent 
hopping |p|-pj|. In addition to the above peculiarities, 
the band/local duality of spin in the DE system further 
results in the following important point: Going to mo- 
mentum representation in Eq. (^) one immediately real- 
izes that the magnon-fermion interaction vertex does not 
vanish at zero-momentum transfer. What is wrong? The 
boson bi should not be regarded as a physical magnon. 
The true Goldstone particle of the DE model is the fol- 
lowing object with both local spin and itinerant features: 
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The composite character of the physical magnon Bi is 
the price one has to pay for spin-charge separation and 
for the rearrangement of the original Hilbert space. The 
itinerant component of Bi is of order 1/S only. However, 
the spin stiffness is itself of the same order, and the 1/S 
correction in Eq. is of crucial importance to ensure 
spin dynamics consistent with the Goldstone theorem. 
Now, commuting Eq. (Q) with Eq. (Q) one finds the mean- 
field magnon dispersion uj p and the correct momentum 
structure of the magnon-fermion scattering vertex: 

[B p , H] = u p B p + ±Y, A f( k ) c L c k- q ,pB p+q , (5) 
q 



Here, uj p = zD s (l - j p ), D s = JbeS, Jde is de- 
fined above, z = 6, and the form factors 7& = 
z^Es^MikS), if - (zty^tfexpiikS). The 
physics behind the second term in Eq. (|^) is the tem- 
poral fluctuations of the "exchange constant" due to the 
charge and orbital dynamics. From now on, the conven- 
tional diagrammatic method with bare magnon-fermion 
vertex given in Eq. (||) can be used. 

Now we turn to the correlations in the fermionic band. 
We separate the charge and orbital degrees of freedom 
by the parameterization c, La — h\fia, where the slave 
boson hi and the fermion /j Q represent the density and 
orbital fluctuations, respectively Further, approxi- 
mating 8(cLcjp) = {h]hi) ft a fjt} + {f}jjf))h]hi, one finds 
from Eq. (H) two different contributions to the magnon 
self energy. 

Superexchange. — At low dopings, the virtual charge 
transfer across the Hubbard gap becomes of importance. 
In limit of large Jh , transitions to the high-spin interme- 
diate state with energy Ui = (V — Jh) <C U', U domi- 
nate. The corresponding superexchange Hamiltonian is 
then obtained to be: 
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Here, J SE = {t 2 /U 1 )[2S(2S + l)]" 1 and t x ' v = 
{u z ± \/3a x ) /2, r z = a z with Pauli matrices cr 2 ^ 2 acting 
on orbital subspace f a . Eq. is the generalization of 
orbitally degenerate superexchange models (|Jl^] for ar- 
bitrary values of spin. Superexchange is of ferromagnetic 
nature because the Jh — * oo limit is assumed, and the 
exchange strength depends on orbital orientations Jl3| . 
In an orbitally disordered state one can represent the 
orbital part of Eq. (||) in terms of fluctuating bond op- 
erators fi a fjf3 |§. Then the effect of superexchange is 
simply to add a new term 
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to the scattering amplitude in Eq. (g|) and to renormalize 
D s — > {Jde + Jse)S — txo{x + x )/2S. Here the constant 
xq = 2xot/U\ with mean-field parameter xo = (/ t u_ c /i) ~ 
h(l — x) defines the effective doping level below which the 
superexchange contribution becomes of importance. 

To summarize up to now: We have developed a quan- 
tum theory of the double-exchange model at T < Tc 
which also includes the correlation effects present in tran- 
sition metal oxides. The obvious advantage of this ap- 
proach over earlier treatments (^^] is that the theory 
has a transparent structure adjusted to describe the low- 
energy spin dynamics, and the large atomic scale Jh does 
not appear in the calculations, either. Yet the corrections 
t/ Jh can be accounted for perturbatively, if necessary. 

Indirect magnon-phonon coupling via orbital sector. — 
An important piece of physics relevant to the experiment 
jjj is still missed in model (|l]), that is the Jahn- Teller 
orbital-lattice coupling This can be written as: 



Hjt = -go [Quel 



hid; 



(8) 



where Q2 and Q3 are the distortions of appropriate sym- 
metry ||. The deformation energy including intersite 
correlations is given by 



with Q* /v = (Q 3 ± V3Q 2 )/2, Q§ = Q 3 , and Q = 
(Q2,Q3)- In general, the JT interaction strongly mixes 
the orbital and lattice dynamics, leading to a fluctuation 
of exchange bonds at low (phonon) frequencies. Here, 
we treat this problem only perturbatively as shown in 
Fig. ||. We arrive at the following effective spin-phonon 
Hamiltonian: 



(9) 



H s -ph - - £ (g Pq Q- q ) B P B p+q . 



(10) 



The coupling constants gjgf = #oGrj(»7g -ri^ —T}j^ q ) / S ', 
where ii q 2) = -\fi{c x - c y )/2, r) q ^ = c z - \c x - 





FIG. 2. Effective spin-phonon coupling. Perturbatively, it 
is controlled by the orbital susceptibility (fermion bubble), 
Jahn- Teller (open dot oc go) and double-exchange (filled dot 
oc t) interaction vertices. 
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FIG. 3. Magnon self energies. 
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j;c y , c a = cos<7 Q , and the parameter ao 

*0)<(// + c/i)K))u,=0- 

Magnon self energies — We are now in the position to 
calculate the renormalization of the magnon dispersion. 
The leading 1/S 2 corrections are shown in Fig. ||. Charge 
and orbital susceptibilities in Figs. [§](a) and ||(b) , respec- 
tively, are calculated using mean-field Green's functions 
in slave boson h and fermion / subspaces. For the spec- 
tral density of JT phonons in Fig. |||(c) we use 



P± h (w,<?) = - — 



7T W ± {Uj-Ujff+T 2 



(11) 



accounting phenomenologically for the damping T due to 
coupling to orbital fluctuations. The phonon dispersions 
w± = A lq ±( K J q + K 2 q )^/ 2 withc ph = {K/M) 1 ' 2 



, . ^1q T ^3q) J ' w"-" w Q 

follow from Eq. (^|). Here, Ki q = 1 + k\(c x 

1 (2) 
K2q = klVq 



K 3q = ki7]q , and k\ — K\jK. The ex- 
pressions obtained from the diagrams in Fig. ^ contain 
summations over momentum space which we perform nu- 
merically. We find the effect of charge fluctuations on the 
magnon spectrum to be quite featureless and moderate 
(see Fig. ||), which is due to the fact that the spectral 
density of charge fluctuations lies well above the magnon 
band. On contrary, the relatively low-energy fluctuations 
of the orbital and lattice degrees of freedom are found to 
affect the spin-wave dispersion in a peculiar way, partic- 
ularly in (0, 0, q) and (0, q, q) directions. Fig. |^ shows the 
magnon dispersion. The hopping amplitude t — 0.4 eV 
is chosen to fit the spin stiffness in Pro.63Sro.37Mn03; 
further we use U\ = 4 eV [ fI2| . The phonon contri- 
bution depends on the quantity (30^0) 2 /2K = Ej^a 2 ,. 
We set E JT al = 0.004 eV [[§, w£ h = 0.08 eV [[TJ, 
r = 0.04 eV. Our key observation is the crucial effect 
of intersite correlations of JT distortions, controlled by 
k\, on the magnon dispersion (see Fig. |^). To explain the 
experimental data Q we are forced to assume ferro-type 
correlations (ki < 0). We interpret this surprising result 
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FIG. 4. Magnon dispersion along (0, 0,£), (£,£,0), and 
(£j£;£) directions, where £ = 0.5 at the zone boundary. Solid 
lines represent the theoretical result including charge, orbital, 
and lattice effects and are fitted to experimental data |l| de- 
noted by circles. For comparison the bare dispersion and the 
one including only charge effects are indicated by long-dashed 
and dashed lines, respectively. The parameter ki — —0.33. 



in the following way: Conventionally k\ > correspond- 
ing to AF order of JT distortions and orbitals is expected 
in undoped compounds |6j. However, in the doped case 
charge mobility prefers a ferro-type local orientation of 
orbitals which minimizes the kinetic energy. This com- 
petition between JT and kinetic energies can be simu- 
lated by tuning k±. At large enough doping, ferro-type 
orbital correlations are expected to prevail, hence effec- 
tively ki < 0. In fact, a ferro-type ordering of orbitals 
leading to a layered AFM spin structure is experimentally 
at dopings about x 



observed 17 18|] at dopings about x — 0.5. As this insta- 
bility is approached, low-energy, quasi static fluctuations 
of exchange bonds develop, yielding a magnon evolution 
as shown in Fig. ||. Remarkably, the small q spin stiffness 
is not affected by this physics, while a strong reduction 
of Tc by soft magnons at the zone boundary is predicted. 
This explains the origin of the anomalous enhancement 
of the D s /Tc ratio in low-Tc manganites |19|. 

In summary, we have presented a theory of spin dy- 
namics in a model relevant to manganites, emphasiz- 
ing particularly the interplay between double-exchange 
physics and orbital-lattice dynamics. The unusual 
magnon dispersion in lower- Tc compounds is explained 
as due to the proximity to orbital-lattice ordering. Ap- 
parently, strongly correlated orbital fluctuations play a 
crucial role in the physics of manganites. 

We would like to thank H. Y. Hwang for stimulating 
discussions. 



100 



80 



> 
E 



60 



k, 

+0.3 
0.0 
-0.3 
-0.33 
-0.333 





FIG. 
Ejtclo 
ki - 



0.25 0.5 0.25 0.5 0.25 
(0,04) (5,4,0) (UA) 

5. Magnon dispersion for different values of fci. 
= 0.006 eV is used. The softening enhances as 
3 corresponding to ferro-type orbital-lattice order. 
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